Measurement of the pressure dependence of 
air fluorescence emission induced by electrons 



AIRFLY Collaboration: M. Ave^, M. Bohacova^, 
B. Buonomo*^, N. Busca'^, L. Cazon^, S.D. Chemerisov'^, 
M.E. Conde^\ R.A. Crowell^ P. Di Carlo^ C. Di Giulio^ 
M. Doubrava^, A. Esposito'^, P. Facal^, F.J. Franchini^, 
J. HorandeP, M. Hrabovsky^, M. larlori*^, T.E. Kasprzyk^, 
B. Keilhauer\ H. Klages-i, M. Kleifges'^, S. Kuhlmann^, 
G. Mazzitelli^ L. Nozka^, A. Obermeier\ M. Palatka^, 

S. Petrera°, P. Privitera^'*, J. Ridky^, V. Rizi^ 
G. Rodriguez^, F. Salamida*^, P. Schovanek^, H. Spinka*^, 
E. Strazzeri^ A. Ulrich^ Z.M. Yusof^, V. Vaceks, 
P. Valente™, V. Verzi^, T. Waldenmaier -i 

University of Chicago, Enrico Fermi Institute, 5640 S. Ellis Ave., Chicago, IL 

60637, United States 

^Institute of Physics of the Academy of Sciences of the Czech Republic, Na 
Slovance 2, CZ-182 21 Praha 8, Czech Republic 

^Laboratori Nazionali di Frascati dell'INFN, INFN, Sezione di Frascati, Via 
Enrico Fermi 40, Frascati, Rome 00044, Ito-ly 

'^Argonne National Laboratory, Argonne, IL 60439 United States 

^ Dipartimento di Fisica deU'Universitd de I'Aquila and INFN, Via Vetoio, 1-67010 

Coppito, Aquila, Italy 

^ Dipartimento di Fisica deU'Universitd di Roma Tor Vergata and Sezione INFN, 
Via della Ricerca Scientifica, 1-00133 Roma, Italy 

^ Czech Technical University, Technicka 4, 16607 Praha 6, Czech Republik 

^ Departamento de Fisica de Particulas, Campus Sur, Universidad, E-15782 

Santiago de Compostela, Spain 

' Universitdt Karlsruhe (TH), Institut fiir Experimentelle Kernphysik (lEKP), 
Postfach 6980, D - 76128 Karlsruhe, Cermany 

^ Forschungszentrum Karlsruhe, Institut fiir Kernphysik, Postfach 3640, D - 76021 

Karlsruhe, Cermany 

^ Forschungszentrum Karlsruhe, Institut fiir Prozessdatenverarbeitung und 
Elektronik, Postfach 3640, D - 76021 Karlsruhe, Cermany 

^Physik Department E12, Technische Universitdt Muenchen, James Franck Str. 1, 

D-85748 Garching, Cermany 

""Sezione INFN di Roma 1, Pie. A. Moro 2, L00185 Roma, Italy 



Preprint submitted to Elsevier 



5 February 2008 



Abstract 



The fluorescence detection of ultra high energy (> 10^^ eV) cosmic rays requires a 
detailed knowledge of the fluorescence light emission from nitrogen molecules, which 
are excited by the cosmic ray shower particles along their path in the atmosphere. 
We have made a precise measurement of the fluorescence light spectrum excited 
by MeV electrons in dry air. We measured the relative intensities of 34 fluores- 
cence bands in the wavelength range from 284 to 429 nm with a high resolution 
spectrograph. The pressure dependence of the fluorescence spectrum was also mea- 
sured from a few hPa up to atmospheric pressure. Relative intensities and collisional 
quenching reference pressures for bands due to transitions from a common upper 
level were found in agreement with theoretical expectations. The presence of argon 
in air was found to have a negligible effect on the fluorescence yield. We estimated 
that the systematic uncertainty on the cosmic ray shower energy due to the pressure 
dependence of the fluorescence spectrum is reduced to a level of 1% by the AIRFLY 
results presented in this paper. 

Key words: Air Fluorescence Detection, Ultra High Energy Cosmic Rays, 
Nitrogen Collisional Quenching 

PACS: , 96.50.S-, 96.50.sb, 96.50.sd, 32.50.-Fd, 33.50.-j, 34.50.Fa, 34.50.Gb 



1 Introduction 



The detection of ultra high energy (> 10^^ eV) cosmic rays using nitrogen flu- 
orescence light emission from extensive air showers (EAS) is a well established 
technique, used by the Fly's Eye P, HiRes [2], and Pierre Auger Observatory 
[3] experiments, and planned for the Telescope Array |4J, which is presently 
under construction. It has also been proposed for the satellite-based EUSO 
[5] and OWL [6] projects. Excitation of atmospheric nitrogen by EAS charged 
particles induces fluorescence emission, mostly in the wavelength range be- 
tween 300 to 430 nm. Information on the longitudinal EAS development can 
be obtained by fluorescence telescopes by recording the light intensity as a 
function of time and incoming direction; this information is related to the pri- 
mary cosmic ray energy and type. However, the fluorescence light yield from 
EAS charged particles must be well known at each point within the shower, 
and corrections applied for atmospheric effects between the shower and the 
telescope for an accurate primary energy determination. Thus, the intensities 
of the fluorescence bands should be measured over a range of air pressure and 
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temperature corresponding to altitudes up to several tens of km, the typical 
elevation of EAS development in the atmosphere. 

A number of other experiments have made measurements of the fluorescence 
light yield pertinent to EAS. These include early low-energy stopped-particle 
results in air by Bunner [7] and electrons in air by Davidson and O'Neil [8]. 
More recently, Kakimoto et al. [9] obtained the light yields with three narrow 
and one broad band optical filters with electrons from a radioactive source at 
1.4 MeV and from a synchrotron at 300, 650, and 1000 MeV in dry air and 
pure nitrogen. Nagano et al. (TU] measured the light yields in pure nitrogen and 
dry air with electrons of average energy 0.85 MeV from a radioactive source, 
using 14 narrow-band filters. Belz et al. [11] used 28.5 GeV electrons in pure 
nitrogen and dry air to measure the light yield through a HiRes broad band 
optical filter (~300-400 nm) as a function of pressure. Colin et al. P3| mea- 
sured the fluorescence light yield with a broad band optical filter at electron 
energies of 1.5 MeV, 20 GeV and 50 GeV. The application of these fluores- 
cence measurements to EAS experiments has been described recently by de 
Souza et al. [13], Keilhauer et al. [11], and Arqueros et al. [T5] . 

The uncertainty on the fluorescence light yield is currently one of the main 
systematic uncertainties on the cosmic ray energy determination by EAS ex- 
periments which employ the fluorescence technique. The absolute fluorescence 
yield is known only at the level of 15% and for a few electron energies. Recent 
spectral measurements are limited by the use of optical filters, while the early 
data of Bunner are based on measurements with coarse spectral resolution 
and large systematic uncertainties. The data reported here are part of a pro- 
gram by the AIRFLY (AIR FLuorescence Yield) collaboration to measure the 
fluorescence light yield with significantly improved precision over the electron 
kinetic energy range from keV to GeV using several accelerators. This paper 
describes the measurement of the pressure dependence of the relative yields 
of the fluorescence bands in the wavelength range 284 - 429 nm. The absolute 
yield, as well as the temperature and humidity dependence of the fluorescence 
spectrum, are currently being analysed, and will be reported elsewhere. 

The relative intensities of 34 fluorescence bands over the wavelength range 
284 - 429 nm were measured in dry air at 800 hPa with a high resolution 
spectrograph. The pressure dependence of the fluorescence spectrum was also 
studied from a few hPa up to atmospheric pressure. The high beam current 
needed for a measurement of the fluorescence spectrum was provided by a 
DC beam of 3 MeV electrons. The pressure dependence of the 337 nm band 
in dry air was measured in a different set-up using a narrow band optical 
filter and a photomultiplier tube. We used a 14 MeV electron beam for this 
measurement, which had the advantage of much better stability of the beam 
position, a small spot size and reduced multiple scattering effects compared 
to the 3 MeV beam. These beam characteristics were important to reduce 
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the systematic uncertainties. We measured the colhsional quenching reference 
pressure of the 337 nm band, p'„j^(337), by studying the ratio of fluorescence 
emission in nitrogen and air. With this method, we ehminated the bias from 
undetected hght due to secondary electrons escaping the detector's field of 
view at low pressures. The colhsional quenching reference pressures of the 
other fluorescence bands were obtained from the pressure dependence of their 
relative intensities, using (337) as normalization. 

This paper is organized as follows. Section 2 presents equations used to param- 
eterize the pressure dependence of the fluorescence light yield. The experimen- 
tal method and hardware are given in Section 3. The measured fluorescence 
spectrum at 800 hPa in dry air, the method used to estimate the band in- 
tensities, and a comparison to theoretically expected intensities are discussed 
in Section 4. The 337 nm pressure dependence measurements at 14 MeV are 
given in Section 5, as well as a new analysis of the data which minimize the 
systematic uncertainties due to the spatial distribution of the fluorescence 
emission induced by secondary electrons. The measurement of the pressure 
dependence of the remaining spectral bands is presented in Section 6. Section 
7 presents applications to EAS, and Section 8 summarizes this work. 



2 Fluorescence yield 



Electrons passing through air excite the nitrogen molecules, directly as well as 
through secondary electrons produced along the path in the gas. Nitrogen de- 
excitation results in a fluorescence spectrum which, in the range 300 to 430 nm, 
consists mainly of transitions from the so-called second positive system (2P) of 
molecular nitrogen N2 and the first negative system (IN) of ionised nitrogen 
molecules N^. Throughout the paper, the notations 2P{v,v') and lN(f,f'), 
corresponding to C^Uuiv) B'^Ug{v') and B^J:+{v) X^J:+{v') [16|, will 
be used to indentify the relevant electronic- vibrational transitions. Not all the 
excited nitrogen molecules emit fluorescence photons, since they may trans- 
fer their energy to other molecules through collision. This quenching process 
introduces a dependence of the fluorescence emission on the gas pressure and 
temperature, because the colhsional rate depends on the average separation 
distance and velocity of the molecules. 

The process of fluorescence emission induced by electrons in nitrogen gas can 
be described in terms of fluorescence efficiency [7] [10], namely the ratio of 
the energy emitted by the excited gas in fluorescence photons to the energy 
deposited in the gas by the electron. The fluorescence efficiency of photons of 
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wavelength A at a given nitrogen pressure p is usually parameterized as: 

*N2(A,p) 



1 + 



where p'-^^{X) is the collisional quenching reference pressure, and (A) is the 
fluorescence efficiency in absence of collisional quenching (in the limit j9 — 
the distance between nitrogen molecules becomes very large and they cannot 
de-excite by collisions). Here, the fluorescence yield is defined as the ratio of 
the fluorescence efficiency to the photon energy E\, that is the number of 
photons emitted by the excited gas per energy deposited by the electron: 

>^N.(A,p) = -^^, (2) 



where FiO.(A) = <f?,,(A)/i?A. 

The quenching reference pressure for nitrogen can be written in terms of the 
lifetime of the excited state to decay to any lower state Tq and the cross section 
for nitrogen-nitrogen collisional de-excitation (Tnn [Z] [ID] : 

where Mn is the N2 molecular mass, k is Boltzmann's constant, and T is 
temperature. 

The fluorescence yield in air is given by: 

>^».(A,P) = ^|^, (4) 
p' (A) 



where /nj is the fraction of nitrogen molecules in air (79%). The quenching 
reference pressure in air takes into account also nitrogen-oxygen collisional 
de-excitation: 



/N20"NN + J02'^^0\\ ^TTT - —, H 3—, \S>) 



Pair \/T^M^kT \ V 2Mo j p'o2 

where /02 is the fraction of oxygen molecules in air (21%), Mq is the O2 
molecular mass and uno is the cross section for nitrogen-oxygen collisional 
de-excitation. In this paper, we will not interpret the measurements in terms 
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of collisional cross sections, but rather use the quenching reference pressure 
for each band as a phenomenological description of the data. Note that the 
fluorescence yield process finds in Eqs. (jl]) - ([S]) its simplest description. Mech- 
anisms which may modify Eqs. dl])-® have been discussed in the literature 
(see for ex. [H] [S] ) . On the other hand, it will be shown in the following 
Sections that the ansatz described by Eqs. dl])-© does indeed provide a good 
description of our data. 

The typical arrangement in fluorescence experiments involves a photon detec- 
tor, e.g. a photomultiplier, viewing a portion of gas volume (N2 or air) along 
the electron path. The detected signal is thus proportional to the number of 
fluorescence photons A^|"* emitted in the detector field of view: 



where E^^^^ is the energy deposited by the electron in the gas volume viewed 
by the photon detector. The energy £"1"^ depends on the specific geometry of 
the experiment, and can be estimated by a detailed Monte Carlo simulation, 
including effects like multiple scattering and secondary electrons. Nevertheless, 
for a given experimental arrangement, we can write: 



where {dE/dX)gas is the collisional energy loss, Pgas is the gas density and 
Dgas{p) is an effective length which takes into account the specific experimental 
arrangement. Notice that a pressure dependence has been included in Dgasip)- 
In fact, as the gas pressure goes down, an increasing fraction of secondary 
electrons escapes the detector field of view, and correspondingly part of the 
fiuorescence emission is not detected. The net effect is a reduction of the 
effective length Dgas, which is thus expected to be a decreasing function of p. 

The number of detected fiuorescence photons can then be written: 



where Rgas is the specific gas constant, T is the gas temperature and the gas 
equation of state has been used. 



(6) 




(7) 
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3 Experimental method 



Measurements were performed at two Argonne National Laboratory accelera- 
tors. Spectra from 284 - 429 nm were recorded with 3 MeV electrons in various 
gases and pressures using the Chemistry Division's Van de Graaff (VdG) elec- 
tron accelerator [19] [20]. The relative light yield of the 337 nm band was 
measured at 14 MeV at the Argonne Wakefield Accelerator (AWA) Facility 
[2T] in dry air and pure nitrogen as a function of pressure. This relative light 
yield, combined with the measured spectra, allowed the pressure dependence 
of the fluorescence bands to be derived. 



3. 1 Van de Graaff Measurements 



The Van de Graaff accelerator is capable of accelerating electrons and protons 
to kinetic energies from 0.5 to 3 MeV, with beam currents in excess of 10 ^A. 
The VdG can operate in DC current or pulsed mode. For the spectrum mea- 
surements it was operated in the DC current mode with typical beam currents 
of ~ 10 /iA, and nominal beam kinetic energy of 3.0 MeV. After exiting the 
VdG, the electron beam was bent 30° in an electromagnet and was focused 
near the exit window 2.00 m from the magnet center. The electrons left the ac- 
celerator vacuum through a 35 mm diameter, 0.152 mm thick dura-aluminum 
window; see Fig. [H An additional quartz window was placed downstream of 
the aluminum window during beam tuning. The beam spot image from the 
quartz was reflected with a mirror to a camera. The beam spot size was typi- 
cally 6 mm diameter, and a side-to-side beam motion of approximately 5 mm 
was observed due to small (< 1%) variations in the VdG energy on time scales 
of seconds. The quartz window and mirror were removed during measurements 
of the spectra. 

The layout for the Van de Graaff measurements is shown in Fig. [H The elec- 
trons entered the pressure chamber after 6 cm of air. Light produced in the 
gas was focused by an aluminum spherical mirror onto the end of a 10 m long, 
1.5 mm diameter pure silica core optical fiber, which brought the light to a 
spectrograph. The optical fiber was placed outside the pressure chamber, and 
light reached the fiber end passing through a quartz window. The spectrograph 
was located behind a concrete block wall and had additional lead shielding to 
protect it from the radiation produced by the VdG. 
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Fig. 1. Layout of the measurements at the Van de Graaff, including the exit win- 
dow from the accelerator vacuum, the pressure chamber, spherical mirror, quartz 
window, and optical fiber. 

3.2 Argonne Wakefield Accelerator Measurements 



The Argonne Wakefield Accelerator can accelerate electrons to kinetic ener- 
gies from 3 MeV to 16 MeV. It was operated in the pulsed mode at 5 Hz, with 
bunches of maximum charge of ~ 1 nC and length ~ 15 ps (FWHM). A new 
set of three quadrupole magnets and two Vernier steering dipole magnets was 
added for this experiment to focus and steer the beam. The nominal beam ki- 
netic energy was 14.0 MeV, with an estimated energy spread of ±0.3 MeV. The 
electrons exited the accelerator vacuum through a 32 mm diameter, 0.13 mm 
thick beryllium window; see Fig. [2J During beam tuning, the beam spot image 
from a quartz plate placed at 45° to the beam direction near the exit window 
was viewed with a camera. The beam spot size was typically 5 mm diame- 
ter, with negligible beam motion. The quartz plate was removed during the 
fluorescence measurements. 

The beam intensity was monitored with an integrating current transformer 
(ICT,[22]), immediately before the beam exit flange. The signal from the ICT 
was integrated, digitized, and recorded for each beam bunch. The ICT was 
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Fig. 2. Layout of the measurements at the Argonne Wakefield Accelerator, including 
the exit window from the accelerator vacuum, ICT, the pressure chamber, quartz 
window, mirror, 337 nm filter, mechanical shutter, and photomultiplier. The last two 
quadrupoles are also shown, but not the third quadrupole or the adjacent dipole. 
Considerable lead shielding of the photomultiplier and mirror were also omitted for 
clarity. 



calibrated by the manufacturer to give an output charge 0.025 times the input 
charge. 



The layout for the AWA measurements is shown in Fig. [21 including the last 
two quadrupoles and ICT. Light produced in the gas propagated out of the 
chamber through a quartz window, mechanical shutter, and 337 nm filter, and 
was reflected by a mirror to a photomultiplier. The shutter could be closed 
remotely to allow measurements of background. A Hamamatsu H7195P pho- 
tomultiplier [23J was used. A second photomultiplier of the same type was used 
to monitor backgrounds. Both were surrounded by considerable lead shielding 
to reduce beam-related backgrounds. The first photomultiplier, mirror, filter, 
and shutter are also shown in Fig. [2l 



A VME data acquisition system was used. The signals from the photomulti- 
pliers and the ICT were integrated by a Lecroy model 1182 charge integrating 
ADC. The accelerator timing signal was used to produce the integrating gate 
of 200 ns width. Signals were recorded for each electron bunch passage in the 
pressure chamber. 
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3.3 Pressure Chamber 



The pressure chamber was constructed of an aluminum tube with various 
flanges welded to it for windows, gauges, gas inlet, and pump-out. The alu- 
minum tube had an inner diameter of 201 mm, length 378 mm, and wall 
thickness 3 mm. The exit window of 0.1 mm thick aluminum was bolted to 
the one end. The entrance window of 0.50 mm thick beryllium, 35 x 55 mm, 
was diffusion bonded to a Confiat flange. All vacuum seals were made with 
0-rings. A top view of the pressure chamber is shown in Fig. [TJ 

A remotely-controlled gas handling and vacuum system was used with the 
pressure chamber. A dry scroll vacuum pump [2lj was used to reduce the pres- 
sure and evacuate the chamber. The chamber pressure was measured [25] at 
the pump-out port. Three types of high purity dry gases were used: 99.9995% 
pure nitrogen, a mixture with argon (78.0% nitrogen, 21.0% oxygen, 1.0% 
argon), and a mixture without argon (79% nitrogen, 21% oxygen). 

3.4 Spectrograph 

Air fluorescence spectra were recorded by an Oriel MS257^*^ spectrograph 
|26] . The main characteristics of the instrument include an asymmetrical 
Czerny- Turner design, F number equal to 3.9, and focal length of 220 mm. 
A holographic grating with 1200 lines/mm and blaze wavelength of 250 nm 
was used. The light from the grating was collected by a 1024x255 CCD pixel 
array (Andor DV420 BU2 [27]). The CCD camera, of 26x26 /im pixel size, 
was back-illuminated and had a large quantum efficiency (~50%) in the 300 
to 400 nm wavelength range. The linear dispersion of spectrograph, with the 
grating used, is 3.2 nm/mm, that gives a wavelength window of about 85 nm. 
This means that the air fluorescence spectrum, which extends over more than 
100 nm, did not fit in a single wavelength window of the spectrograph. Thus, 
in a fluorescence run the spectrograph collected data in sequence, first in the 
range 284-369 nm and then in the range 344-429 nm. For each wavelength 
range, 50 spectra of 1 second integration time were taken. The spectrograph 
operation and data acquisition was fully automated. 

The spectral bands emitted by a mercury pencil lamp (Oriel no. 60635 [26j) 
in the region 290 to 410 nm were used to calibrate the spectral response of 
the spectrograph. The calibration was then refined by including the positions 
of the nitrogen emission bands as seen in the air fluorescence spectrum mea- 
surement. 

A quartz tungsten halogen lamp (Oriel no. 63350 [2^) was used to calibrate 
the relative spectrograph sensitivity as a function of wavelength. The spectral 
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irradiance of the lamp is traceable to NIST [28], and its uncertainty amounts to 
2.5% in the wavelength range of interest for the air fluorescence measurement. 

The spectrograph calibration was performed on the beam line, between air 
fluorescence measurements. A flange on the top of the pressure chamber al- 
lowed the insertion of the calibration lamp along the beam path. In order to 
take into account the spatial distribution of the air fluorescence emission in 
the experiment, the calibration lamp was placed in different positions along 
the beam axis and the recorded lamp spectra were appropriately combined. 
By comparison with the known lamp spectrum, calibration factors for the 
relative intensity as a function of wavelength were determined. The two spec- 
trograph ranges used for the air fluorescence measurement, 284-369 nm and 
344-429 nm, were independently calibrated. Their relative normalization, af- 
ter the calibration procedure, was checked by comparing the intensities of the 
mercury pencil lamp lines emitted in the overlap region of the two spectro- 
graph's ranges. A 2% difference in the measured intensities was found, which 
was taken as an estimate of the systematic uncertainty on the spectrograph 
relative calibration procedure. 



4 Measurement of air fluorescence spectrum 

The fluorescence spectrum in air-like mixtures has been measured at the Van 
de Graaff . Details of the beam and apparatus were given in Section 13. 1[ Mea- 
surements were performed at a gas pressure of 800 hPa and temperature of 
293 K. A blow-up of the spectrum in the wavelength region between 343 nm 
and 360 nm is shown in Fig. [31 Notice that the observed lineshapes show an 
asymmetric tail towards smaller wavelengths, due to the convolution of the ro- 
tational levels with the finite resolution of the spectrograph. Thus, in general, 
the integral counts associated with a given band can have contributions also 
from nearby bands. Rather than attempting to disentangle individual band 
integrals, which would result in additional systematic uncertainties, the band 
intensities were estimated by integrating the CCD counts in a wavelength in- 
terval around each band. An example of the procedure is shown in Fig. [3J 
Wavelength intervals up to several nm were used, depending on the distance 
of nearby bands and the amount of rotational contribution. In any case, most 
of the signal of each band is contained in only 2 nm. In the following, the 
wavelength of the main band is used to identify the corresponding wavelength 
interval. The measurement of the light intensity in wavelength intervals of nm 
size is well suited for the purpose of calibration of the cosmic ray fluorescence 
detection technique, and assigning the full intensity to the main band of the 
wavelength interval will result in a negligible systematic uncertainty. 

The measured fluorescence spectrum contained a small level of background. 
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Fig. 3. Example of band integrals. Horizontal lines indicate the wavelength ranges 
used for the integration of the band fluorescence signal. The main fluorescence band 
corresponding to each interval is also given. 

which was estimated in selected wavelength regions where fluorescence bands 
were not present. For each of these background regions, the average number of 
background counts per wavelength bin was determined. In order to estimate 
the background level under the fluorescence bands, a linear interpolation was 
performed between the measurements of consecutive background wavelength 
regions. The procedure gave an estimate of the average number of background 
counts in each wavelength bin over the full wavelength region. The background 
level was found to be approximately constant. The amount of background un- 
der each fluorescence band was determined by summing the estimated average 
number of background counts for each bin of the wavelength interval associ- 
ated to the band. The band intensity was then determined by subtracting the 
corresponding background. 

Several fluorescence bands (346.3, 350.0, 353.7, 357.7, 366.1, and 367.2 nm) 

were common to the 284-369 nm and the 344-429 nm spectrograph's ranges. 
The measured intensity of these bands was used to normalize the two spectro- 
graph's ranges. In fact, any change in the beam intensity or position during 
the short time (30 s) needed by the spectrograph to move from the first to the 
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second wavelength range would result in a difference in the measured intensity 
of the common bands. The intensity of the common bands was found to be 
2% lower in the 344-429 nm range, and the intensity of all measured bands 
in this range were corrected correspondingly. The smallness of the correction 
illustrates the good stability of the beam and the reliability of the relative 
calibration of the spectrum. 



The measured spectrum with the 78% N2 - 21% O2 - 1% Ar mixture is shown 
in Fig. m The measured intensities, Ja, relative to the intensity of the 337 nm 
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A (nm) 


A interval (nm) 


h (%) 


296.2 


292.5-297.0 


5.16 ±0.29 


366.1 


363.2-366.4 


1.13 ±0.08 


297.7 


297.0-299.6 


2.77 ±0.13 


367.2 


366.4-367.5 


0.54 ± 0.04 


302.0 


301.5-303.3 


0.41 ±0.06 


371.1 


367.6-371.7 


4.97 ±0.22 


308.0 


306.8-309.3 


1.44 ±0.10 


375.6 


371.7-376.3 


17.87 ±0.63 


311.7 


309.3-312.3 


7.24 ±0.27 


380.5 


376.3-381.4 


27.2 ± 1.0 


313.6 


312.3-314.1 


11.05 ±0.41 


385.8 


383.0-386.0 


0.50 ±0.08 


315.9 


314.1-316.7 


39.3 ± 1.4 


387.7 


386.0-388.0 


1.17 ±0.06 


317.6 


317.0-318.4 


0.46 ± 0.06 


388.5 


388.0-388.7 


0.83 ± 0.04 


326.8 


325.6-327.1 


0.80 ±0.08 


391.4 


388.7-392.1 


28.0 ± 1.0 


328.5 


327.1-329.0 


3.80 ±0.14 


394.3 


392.1-394.9 


3.36 ±0.15 


330.9 


329.0-331.3 


2.15 ±0.12 


399.8 


394.9-400.5 


8.38 ±0.29 


333.9 


331.3-334.3 


4.02 ±0.18 


405.0 


400.5-406.6 


8.07 ±0.29 


337.1 


334.3-338.4 


100.00 


414.1 


412.5-414.4 


0.49 ± 0.07 


346.3 


344.2-347.2 


1.74 ±0.11 


420.0 


416.6-420.6 


1.75 ±0.10 


350.0 


347.2-350.6 


2.79 ±0.11 


423.6 


420.7-424.0 


1.04 ±0.11 


353.7 


350.6-354.4 


21.35 ±0.76 


427.0 


424.0-427.4 


7.08 ±0.28 


357.7 


354.4-359.9 


67.4 ±2.4 


427.8 


427.4-428.6 


4.94 ±0.19 



Table 1 

Measured fluorescence band intensities in dry air at 800 hPa pressure and 293 K 
temperature. The intensity of the 337 nm band was used for normalization. The 
wavelength intervals used for the signal integration are also reported. 



band, of 34 fluorescence bands are reported in Table [H, together with the 
corresponding wavelength intervals used for the signal integration. The fluo- 
rescence spectrum was measured ten times, and from the observed r.m.s. of 
the ten measurements of each band intensity, uncertainties were estimated. 
To obtain the total uncertainty quoted in Table [H a systematic uncertainty 
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of 3.5%, assigned to the calibration of the relative spectrograph sensitivity, 
was summed in quadrature. The systematic uncertainty of the relative spec- 
trograph sensitivity accounts for the uncertainty on the absolute calibration 
of the lamp emission (2.5%), the difference in the mercury pencil lamp line 
intensities measured in the two wavelength ranges of the spectrograph (2%), 
and the uncertainty associated with the spatial distribution of the fluorescence 
light source (1.5%). The fluorescence light is in fact emitted along the beam 
axis, and the calibration procedure of relative spectrograph sensitivity (cf. Sec- 
tion 13. 4p derived calibration factors from measurements with the calibration 
lamp placed in different positions along the beam axis. If calibration factors 
obtained only from a measurement with the calibration lamp placed at the 
optics center are used, the relative fluorescence band intensities changed by 
at most 3%. Half of this shift was conservatively taken as an estimate of the 
associated systematic uncertainty. 

Several checks were performed. The linearity of the fluorescence emission with 
beam currents from 0.2 to 14 /iA was verifled. Possible systematic effects due to 
the beam position were investigated. The fluorescence spectrum was measured 
with the beam moved ±1 cm in the directions transverse to the nominal beam 
axis. No difference beyond the statistical uncertainty in the relative intensities 
of the fluorescence bands was found. Several models for background evaluation 
were tested, which always resulted in changes of the relative band intensities 
within the quoted uncertainties. A measurement of the fluorescence spectrum 
in pure nitrogen gas was performed, which showed that all the observed bands 
in the air fluorescence spectrum are associated with nitrogen excitation. 

In order to assess the relevance of argon on air fluorescence, the fluorescence 
spectrum emitted by a 79% N2 - 21% O2 gas mixture was measured at the same 
pressure and temperature conditions as the mixture with argon. Fig. O shows 
the correlation of the relative intensities of the 34 fluorescence bands measured 
with the 79% N2 - 21% O2 and the 78% N2 - 21% O2 - 1% Ar mixture. A linear 
flt yielded an intercept consistent with zero and a slope consistent with unity, 
within a few per mil, which indicates that the two spectra are compatible to 
a high degree of accuracy. On the basis of this observation, the effect of argon 
on air fluorescence at atmospheric pressure is completely negligible. 

The 2P and IN bands of nitrogen emission can be easily identifled in our 
measured fluorescence spectrum. A group of weak bands (302.0, 308.0, 317.6, 
346.3, 366.1, 387.7 nm) is consistent with the Gaydon-Herman (GH) bands 
[IS] [2n]. The GH(0,l-3) bands at 296.7, 311.9 and 328.3 nm are also present, 
with a large overlap with the 2P(3,l-3) bands at 296.2, 311.7 and 328.5 nm. 
The excitation of the Gaydon-Herman bands by electrons in air at atmo- 
spheric pressure has already been observed |30]. Within a vibrational band 
system, the relative intensities are expected to be equal to the ratio of the 
corresponding Einstein coefficients, and can be calculated [21] [32] • However, 
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Fig. 5. Effect of argon on the air fluorescence yield. Relative intensities of the 34 
fluorescence bands measured in a 79% N2 - 21% O2 mixture vs the corresponding 
intensities measured in the 78% N2 - 21% O2 - 1% Ar mixture. A linear fit of the 
data is also shown, which indicates that the yield in the two gas mixtures is the 
same to a high level of accuracy. 



the intensities reported in Table [T] cannot be directly compared with the the- 
oretical expectations, due to the contamination from nearby bands. A simple 
extrapolation of the band's spectral shape was used to estimate the contami- 
nation, and intensities were corrected accordingly. The corresponding relative 
intensities for the 2P(0,i), 2P(l,i), 2P(2,i), and lN(0,i) bands are reported in 
Table together with the theoretical predictions [21] [22] • Only bands where 
the contamination could be reliably estimated are included in Table [2l The 
measured relative intensities are close to the theoretical expectations and this 
represents an independent cross-check of the relative spectrograph sensitivity 
calibration. 
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Table 2 

Measured relative intensities within bands and theoretical expectations based on 
the ratio of Einstein coefficients. 



5 Pressure dependence of the 337 nm band 

A precise measurement of the pressure dependence of the 337 nm band has 
been performed at the AWA facility. The experimental set-up has been de- 
scribed in Section 13. 2[ A narrow band interference filter centered at 337 nm 
was used for this measurement. The filter transmission curve is shown in Fig. 
ini together with the fluorescence emission bands as measured by the spec- 
trograph (Section HI). The filter transmission was measured by three different 
groups of the AIRFLY collaboration, which yielded consistent results at the 
percent level. The filter selects mainly the 337 nm band, the contamination 
from neighbouring bands being only 1.7%. 

In order to check linearity and possible pedestal changes within a run, the 
electron beam was operated in a mode allowing the bunch charge to fluctuate 
over a wide range. An example of correlation of the photomultiplier signal 
SpMT and the ICT signal Sjct is shown in Fig. [71 for a run at 180 hPa in pure 
nitrogen. A linear fit Spmt = SpL ■ Sjct + b to the data was performed, and 
the fitted slope Spp was taken as an estimator of the fluorescence signal. The 
same procedure, applied on data taken with the shutter closed in front of the 
photomultiplier tube, allowed an estimate of the background level. A typical 
background correlation to the beam intensity is also shown in Fig. [71 Several 
background runs were taken during the pressure scans. In addition, during the 
fluorescence runs, the background was monitored by a second photomultiplier. 
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Wavelength (nm) 

Fig. 6. Transmission of the interference filter used to select the 337 nm band (dashed 
line). The measured fluorescence spectrum in arbitrary units is shown as a full line. 



and found to be stable within its statistical uncertainty. The background slope 
was subtracted from the signal slope, and in the following Spi will refer to 
the background subtracted signal. 



Measurements were performed with dry air-like mixtures and pure nitrogen. 
The gas pressure in the fluorescence chamber was varied in the range from 
2 to 1000 hPa, and for each pressure p a value Sfl{p) was measured. The 
gas temperature was 293 K. Only a couple of minutes for each pressure point 
was necessary to collect enough statistics, thus minimizing the possibility of 
systematic changes of the beam and the gas condition, as well as the detector 
response. 

The measured fluorescence signal Sp^Kp) in nitrogen or air is proportional to 
(cf. Eq. (ED): 

S^/lip) ^ Ni- = (§) ^D,Up)^^^^. (9) 

\ / gas gas p'gasi^^^) 
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Fig. 7. Correlation between the photomultiplier signal and the beam ICT signal: 
a) fluorescence run in pure nitrogen at 180 hPa, b) background run. 

where = 1 for pure nitrogen. Notice that the fluorescence emission is 
essentially due to excitation by secondary electrons produced in the gas by 
the beam electrons [33] [M]. As the gas pressure goes down, an increasing 
fraction of secondary electrons does not stop in the photomultiplier fleld of 
view, and correspondingly part of the fluorescence emission is not detected. 
Neglecting the losses due to the secondary electrons escaping the fleld of view 
would cause an overestimation of the quenching reference pressure p' [3T 



In order to avoid this bias, we have studied the ratio of the 337 nm fluorescence 
signal in pure nitrogen to the signal in air, r^^ = Spl{p)/Spl{p). The effective 
length D{p) is expected to cancel in the ratio, since secondary electron inter- 
actions are similar in nitrogen and air thanks to their close molecular masses. 
This was verifled with a detailed GEANT4 [SB] simulation of the experiment. 
The pressure dependence of the fluorescence signal ratio thus becomes: 



(dE\ 



l+p 



+ 



fOo 



p' (337) ' p' (337) 



( — ) 

\d.X J air 



1 + 



(10) 



Pn2(337) 



where /nj and /oj are the fraction of nitrogen and oxygen molecules in the air 
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Fig. 8. Pressure dependence of the ratio of the 337 nm band signal in nitrogen to 
the signal in dry air. The full line is the result of a fit to the full pressure range 
described in the text. 



mixture. The coUisional energy loss (dE/dX) in nitrogen and air at 14 MeV 
was calculated with the method described in [35] . 

The measured pressure dependence of the ratio r^^ of the 337 nm band is 
shown in Fig. [HI The 78% N2 - 21% O2 - 1% Ar air mixture was used for this 
measurement. The full line is the result of the best fit of Eq. fllOl) to the data 
in the range 2-1000 hPa, with p'^^ and p'q^ as free parameters. The best fit 
yielded: 



p;^^(337) = 103.7 ±2.7 hPa, 
p'q^(337) = 3.796 ± 0.076 hPa, 



(11) 
(12) 



where the quoted uncertainties are statistical only, and a x^/n.d.f. = 0.9. 
From Eq. (jSj), a value of p'^^^. for the 78% N2 - 21% O2 - 1% Ar mixture was 
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Fig. 9. Ratio of the fluorescence signal for air-like mixture with 1% Ar to the mixture 
without Ar function of pressure. 

estimated: 



p'^i^{337) = 15.89 ± 0.33 hPa, 



(13) 



where the statistical uncertainty was calculated from the covariance matrix of 
the fitted p'^^ and p'q^ parameters. 

Notice that the contribution of argon to the quenching of the 337 nm band 
was neglected in Eq. ffTOj) . This is justified by the comparison of the pressure 
scan measurements performed with two air-like mixtures. The ratio rAr of the 
fluorescence signal measured with the 78% N2 - 21% O2 - 1% Ar mixture to 
the one with the 79% N2 - 21% O2 mixture as a function of pressure is shown 
in Fig. O A linear fit to the data yielded an intercept of 1.011 ± 0.006 and 
a slope consistent with zero at the level of 10~^, with a x^/n.d.f. = 1.0. The 
effect of argon to the quenching of the 337 nm band is thus negligible at all 
pressures. 

Several checks were performed in order to estimate the systematic uncertain- 
ties on the measured p'. The stability of the photomultiplier gain was moni- 



21 



tored during the measurements by a blue LED. The PMT gain was found to 
be stable at the 1% level, with negligible effect on the determination. The 
range used to fit the fluorescence signal Spmt as a function of the ICT signal 
SicT was varied over a wide interval. The lower limit in the fit of the pressure 
dependence was changed from 2 hPa to 30 hPa. The background was changed 
between its minimum and maximum measured value. The observed shifts in 
the p' value due to these variations were taken as estimates of the system- 
atic uncertainties. A systematic uncertainty of 0.1 hPa was assigned from the 
uncertainty in the measurement of the absolute pressure. The different con- 
tributions to the systematic uncertainty, together with their quadrature sum, 
are summarized in Table [31 The total uncertainty on the measured p', given 
by the quadrature sum of the statistical and systematic uncertainties, is also 
quoted. 



Source 


^P'^, (hPa) 


^P'o, (hPa) 


Ap^, (hPa) 


slope fit range 


3.6 


0.10 


0.50 


background 


2.1 


0.09 


0.40 


pressure fit range 


0.3 


0.01 


0.03 


absolute pressure 


0.1 


0.10 


0.10 


TOTAL SYST. 


4.2 


0.17 


0.65 


TOTAL UNCERTAINTY 


5.0 


0.19 


0.73 



Table 3 

Estimated systematic uncertainties on the quenching reference pressure of the 337 
nm fluorescence band. The total uncertainty, given by the the quadrature sum of 
statistical and systematic uncertainties, is also given in the last line. 

The method of the nitrogen to air ratio has the advantage of eliminating the 
bias due to secondary electrons escaping the field of view of the photomul- 
tiplier. On the other hand, a direct interpretation of the data with Eq. ([9]) 
allows a valuable check of our understanding of the experiment. For this pur- 
pose a detailed GEANT4 [36] simulation of the experiment was performed. 
The PENELOPE [37] package in GEANT4 was used, which allows the track- 
ing of secondary electrons down to 250 eV. In the simulation, the fluorescence 
emission was assumed to be proportional to the energy deposited by the par- 
ticles in the gas. The effective length, which takes into account the field of 
view losses, was found to be well parameterized by D{j)) = DiqoqF^p) from 
the Monte Carlo simulations, where -Diooo is the effective length at 1000 hPa 
and F{p) = (p/lOOO)^ ^^''. The pressure dependence of the fluorescence signal 
can thus be fltted with the function: 

S'fUp) = C ^ , F{p), (14) 
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leaving C and as free parameters. A direct fit of the pressure dependence 
of the 337 nm band fiuorescence signal 5*^1^ in nitrogen with Eq. yielded: 

p' (337) = 101.0 ± 0.6 hPa, (15) 



while a fit of the pressure dependence of the 337 nm fiuorescence signal S' 
in the 78% N2 - 21% O2 - 1% Ar mixture gave: 



air 
FL 



p^i^.{337) = 15.34 ± 0.53 hPa, (16) 

in good agreement with (fTTj) and (fT3|) . To avoid the difficulties of simulations 
at very low gas densities, we limited the fit to a minimum pressure of 15 
hPa. The fact that two different ways of taking into account the effect of 
secondary electrons escaping the field of view, either by cancelation in the 
nitrogen to air fiuorescence ratio or by a parameterization of the light losses 
with simulations, yielded consistent values of p' strengthen the confidence in 
our quenching reference pressure measurements. Also, this consistency shows 
that the simple ansatz of Eqs. dl])-© does indeed provide a good description 
of our data. 

It should be stressed that the effect of secondary electrons escaping the de- 
tector field of view cannot be neglected in the measurement of p'. In fact, a 
fit of the data with F{p) = 1 in Eq. ([HD gave p'^^i337) = 115.5 hPa and 
p^j^(337) = 19.95 hPa, showing a significant bias towards larger p' values. 



6 Pressure dependence of the spectrum bands 



The pressure dependence of the air fiuorescence spectrum bands has been 
measured at the Van de Graaff. The experimental arrangement and beam 
conditions have been described in Sections [3] and HI The air fluorescence spec- 
trum was measured for gas pressures in the range from 4 to 1000 hPa, at 
temperature 293 K. For each pressure, the spectrum band intensities relative 
to the 337 nm band were measured following the procedure described in Sec- 
tion HI From Eq. ([9]), the pressure dependence of the relative intensity of a 
band of wavelength A is given by: 

hip) = C^-^^^. (17) 



(A) 



The measured I\{p) are shown in Figs. [T0lfT2] for one of the ten pressure scans 
which were performed. We could not measure the pressure dependence for 



23 



cs, 0.07 

as 

0.06 



0.05 
0.04 
0.03 




10 



0.04 



2 3 
10 10 

p (hPa) 




p (hPa) 




p (hPa) 



o 

0.025 


1 INI 


Ml 1 1 1 1 M M 1 1 1 1 M M - 






0.02 






0.015 








2 3 
10 10 10 



p (hPa) 



0.03 - 



0.02 



0.01 



0.16 



0.14 - 

0.12 - 

0.1 - 

0.08 — 



1 



0.05 



00 



0.04 



0.03 



0.02 



S 0.05 r 

" 0.045 r 

0.04 r 

0.035 r 

0.03 — 
1 



10 



10 



10 10 

p (hPa) 



10 10 

p (hPa) 




10 10^ 10^ 



p (hPa) 




2 3 
10 10 10 

p (hPa) 



Fig. 10. Pressure dependence of the air fluorescence bands shown in increasing order 
from 296.2 nm to 333.9 nm. Intensities are relative to the 337 nm band. Errors are 
statistical only. The full line is the result of a fit described in the text. 
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Fig. 11. Pressure dependence of the air fluorescence bands shown in increasing order 
from 346.3 nm to 387.7 nm. Intensities are relative to the 337 nm band. Errors are 
statistical only. The full line is the result of a fit described in the text. 
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Fig. 12. Pressure dependence of the air fluorescence bands shown in increasing order 
from 388.5 nm to 427.8 nm. Intensities are relative to the 337 nm band. Errors are 
statistical only. The full line is the result of a fit described in the text. 
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nine of the bands in Table [H due to their low intensity. The full line is the 
result of the best fit of Eq. (fTTl) to the data, leaving C and p'airW as free 
parameters. The 337 nm quenching reference pressure p^j^(337) was fixed to 
the value of 15.89 hPa obtained in Section O (Eq. (IT^ ). In principle, the 
measured /a (800 hPa) given in Table [1] could be used as a constraint for the 
C parameter, therefore reducing the error on p'airi^)- On the other hand, the 
measured values of p'„j^(A) and /a(800 hPa) would become correlated. Leaving 
C free in the fit also minimizes systematic uncertainties, since the calibration 
of the relative spectrograph sensitivity does not affect in this case the measured 
quenching reference pressure. 

Notice that Eq. (ITTj) describes a model where the pressure dependence of a 
spectrum band intensity is given by only one quenching reference pressure. This 
model may not be true in general because more than one band could be 
contained in the wavelength interval used for the signal integration. From 
the quality of the fits, however, it appears that the assumption of a single 
reference pressure is a good approximation. The largest deviation from the 
model, < 10%, is found for the 311.7 nm band, which is expected to have 
contributions of similar size from the GH(0,2) and 2P(3,2) bands. A more 
elaborated model of pressure dependence of the spectrum bands could be used, 
but since it would not change significantly the application of our measurements 
to the fluorescence detection of ultra high energy cosmic rays, the simplest 
approach of Eq. (fTTl) was adopted. 

In order to study possible systematic effects, the measurement was repeated 
several times. Measurements were performed in two separate test beam peri- 
ods, under different beam conditions and light collection optics. Both the 78% 
N2 - 21% O2 - 1% Ar and the 79% N2 - 21% O2 air-like mixtures were used. All 
measurements were found to be compatible within the quoted uncertainties. 

The measured p'airW^ obtained from the weighted average of ten measure- 
ments, are reported in Table HI Two uncertainties are quoted: the first is the 
r.m.s. of the ten measurements, and the second comes from changing the value 
of p^j^(337) in the fit by its total uncertainty Ap^^(337) = ±0.73 hPa (see 
Table ED . 

Spectrum bands within a system should have the same quenching reference 
pressure. This may not be true in our measurements, since for each band the 
signal is integrated over a wavelength interval, which can contain also other 
bands. However, most of the bands are well separated, and contamination is 
in general small. Indeed, we observe that most of the bands of a given system 
have the same p' within errors. Bands which deviated from this behavior, like 
the 387.7 and 427.0 nm bands, are significantly contaminated by nearby bands. 
This good agreement gives further confidence on the quality of our pressure 
dependence measurements. 
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Table 4 

Collisional quenching reference pressures in dry air at 293 K. The quoted uncertain- 
ties are explained in the text. 

7 Application to the fluorescence detection of ultra high energy 
cosmic rays 



The number of fluorescence photons emitted at a given stage of a cosmic ray 
shower development, i.e. at a given altitude h in the atmosphere, is given by: 

iVf = Etr^{h)Ya„{\p.T), (18) 



where is the energy deposited by the shower charged particles in 

the air volume, p and T are the air pressure and temperature at the altitude h. 
Notice that we assumed in Eq. (fT8|) that the fluorescence emission is propor- 
tional to the energy deposited in the gas, which has been tested experimentally 
with good precision [ID] [12] [38] . 

Fluorescence telescopes measure N^°'^'^'^{h) to detect ultra high energy cos- 
mic rays. If the fluorescence yield Yair{X,p,T) is known, the energy deposited 
by the cosmic ray shower can be determined. Integrating along the shower 
path in the atmosphere will give the electromagnetic energy of the shower. 
The fluorescence efficiency must thus be known over a range of pressure and 
temperature corresponding to altitudes up to several tens of km. 
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It is convenient to express the fluorescence efficiency in terms of tlie pliysi- 
cal quantities measured in tliis paper, namely the air fluorescence spectrum 
relative band intensities Ix{po,Tq) (Section H]) and the quenching reference 
pressures p^j^,(A,To) (Sections [MSD- The air pressure po = 800 hPa and tem- 
perature To = 293 K used in the measurements have been introduced in the 
notation for consistency with what will be used in the following. The quench- 
ing reference pressure (cf. Eq. (|5])) can be written in a general way as: 



K.(A,T) = -^-^, (19) 



where H\{T) has been introduced to take into account a possible tempera- 
ture dependence of the coUisional cross sections. In terms of the measured 
quenching reference pressures p^j^(A,To), Eq. (|T9l) becomes: 

pL.(A,T) =pL,(A,To)./^:^^. (20) 

V ^0 J^\\J- ) 



The fluorescence yield can thus be expressed in terms of the measured physical 
quantities: 

y,,,(A,p,T) = y,,,(337,po,ro)JA(po,ro)^ (21) 



This parameterization, together with the values of Tables [H IH and [5l can be 
used by ultra high energy cosmic rays experiments which employ the fluores- 
cence detection technique. In Table |5] our best estimates for the unmeasured 
quenching reference pressure of the nine weak bands not present in Table H] 
are reported. We assumed in Table [5] that all the Gaydon-Herman bands have 
the same p' of the GH(0,4) band. For the unmeasured p' of the 2P(3,i) and 
2P(4,i) bands we took the average value of the measured 2P(3,i=l,3) with a 
large uncertainty. The quenching reference pressure of the 1N(1,2) was taken 
to be equal to the measured p' of the 1N(1,1). 

Notice that in order to fully determine the fluorescence yield, the absolute 
fluorescence yield of the 337 nm band, Yair{'i'iT ,Pq,Tq), and the temperature 
dependence of the coUisional cross sections, H\{T) must be known. Measure- 
ments of these quantities will be the subject of subsequent papers. 

The uncertainties on the relative band intensities and quenching reference 
pressures will result in an uncertainty on the cosmic ray shower energy. For 
a proper estimate of this uncertainty, the spectral sensitivity of the fluores- 
cence detector and the wavelength dependent atmospheric attenuation must 
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Bajiid. 


A (nin) 




Ba,nd. 


A (nin) 




2P(4,4) 


OZD.O 


1 n 1 c 1 1 

ly ± ± 1 


uM(6,2j 


3U2.U 


o 1 1 in 1 1 
zl ± 10 ± 1 




OOO.O 


1 n 1 rr 1 -1 


GH(5,2) 


308.0 


21 ± lOib 1 


2P(3,5) 


367.2 


19 ± 5 ± 1 


GH(6,3) 


317.6 


21 ± lOib 1 


2P(3,7) 


414.1 


19 ± 5 ± 1 


GH(0,5) 


366.1 


21 ± lOib 1 


1N(1,2) 


423.6 


3.9 ± 1.7 ±0.3 





Table 5 

Collisional quenching reference pressures in air at 293 K adopted for the nine weak 
bands for which p' could not be measured. 

be taken into account because they change the relative weights of the fluo- 
rescence spectrum bands. A detailed study of these effects goes beyond the 
scope of this paper. Nevertherless, an estimate may be obtained from the total 
fluorescence yield: 

Ytot{p^T) =Y.Ya„{\p,T), (22) 

A 

where the sum goes over the 34 wavelengths of Table [U In Eq. (1221) we assumed 
H\{T) = 1, i.e. collisional cross sections have no temperature dependence. 
From Eq. f|T8|) . the relative uncertainty on the energy deposited by the cosmic 
ray shower is equal to the relative uncertainty on the total fluorescence yield, 
^Ytotl^tot- In Fig- [131 the relative uncertainty on the total fluorescence yield 
is shown as a function of altitude. The U.S. 1976 Standard Atmosphere [31] 
was used to calculate the pressure and temperature at a given altitude. In the 
calculation of (JYtoti uncertainties on the relative intensities of the spectrum 
bands and the quenching reference pressures were treated as uncorrelated, 
with the exception of the uncertainties related to p^j^(337) which were fully 
correlated for all the reference quenching pressures. Ultra high energy cosmic 
rays deposit most of their energy well below 25 km altitude, which corresponds 
to a vertical slant depth of only 14 g/ cm^. Thus, we may estimate from Fig. [T3] 
that the uncertainties on the measurements reported in this work will result 
in a relative uncertainty on the cosmic ray shower energy at the level of 1%. 

From our knowledge of the spectrograph's calibration procedures and several 
checks performed, we estimate that the systematic uncertainty (3.5%) on the 
relative intensities of the spectrum bands are largely uncorrelated. Even taking 
half of this uncertainty as fully correlated when calculating (JYtot-i relative 
uncertainty on the total fluorescence yield would still be very small, of the 
order of 2%. 

The AIRFLY fluorescence spectrum is significantly more precise than previous 
measurements, which had coarser spectrograph resolution [7] or made use of 



30 




5 10 15 20 25 30 

Altitude (km) 



Fig. 13. Relative uncertainty on the total fluorescence yield as a function of altitude. 



narrow band optical filters [9] fn]\. For the purpose of comparison, the total 
fluorescence yield in the range 300 to 400 nm can be used, since this quantity 
was measured also by [11] and [12] with a broad band optical filter. 

The ratio of the total fluorescence yield in the range 300 to 400 nm as measured 
by other experiments to the one measured by AIRFLY is shown as a function 
of altitude in Fig. [TH All yields were normalized to have the same value at 
ground level. Note that optical atmospheric attenuation and scattering effects 
have not been included in Fig. [Ml Thus, the impact on the primary cosmic 
ray energy cannot be directly derived from this figure. Nevertheless, the total 
yields of Bunner [7] and Belz et al. [H] are in very good agreement with 
AIRFLY. Notice that the measurement of Belz et al. should not suffer from 
systematic uncertainties due to secondary electrons, since they measured the 
pressure dependence of the fluorescence lifetime. Kakimoto et al. p]|, Nagano 
et al. [To] and Colin et al. [12] present significant deviations, which reflect the 
systematically higher values of p' measured by these experiments. 
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Fig. 14. Ratio of tlic total fluorescence yield in the range 300 to 400 nm as measured 
by other experiments to the one measured by AIRFLY as a function of altitude: a) 
Belz et al. b) Bunner c) Nagano et al. d) Kakimoto et al. e) CoUn et al. 

8 Summciry 

We have made a precise measurement of the emission spectrum of nitrogen 
molecules excited by 3 MeV electrons in dry air. We could measure 34 fluo- 
rescence bands in the wavelength range from 284 to 429 nm. The 2P and IN 
systems of molecular nitrogen were found to dominate the fluorescence emis- 
sion. The high resolution spectrograph allowed the identification of a group 
of weaker bands, which was found to be consistent with the Gaydon-Herman 
bands. The relative intensities of bands corresponding to transitions from a 
common upper level were in good agreement with theoretical expectations 
based on the ratio of Einstein coefficients. 

The pressure dependence of the fiuorescence spectrum was also measured from 

a few hPa up to atmospheric pressure. Particular care was taken to avoid the 
bias from undetected light due to secondary electrons escaping the detector's 
field of view at low pressures. For this purpose, the relative intensity of the 
337 nm band was measured as a function of pressure in both nitrogen and air. 
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and from the ratio of intensities measured with the two e; measurement 
of the colhsional quenching reference pressure, p'^j^(337), was obtained. The 
measurement of the fluorescence spectrum at different pressures allowed the 
determination of the colhsional quenching reference pressures of 24 fluores- 
cence bands, in addition to the 337 nm band. Systematic uncertainties were 
minimized by taking the intensities relative to the 337 nm band, and by using 
the value of ^[^^^(337) previously measured. Consistent values of the colhsional 
quenching reference pressures were found for bands due to transitions from a 
common upper level. 

The effect of argon on the fluorescence yield was also investigated, and found 
to be neghgible. 

The high resolution spectra recorded by the spectrograph allowed many more 
closely spaced bands to be separated than in previous experiments. In the 
ratio of fluorescence intensities, which was extensively used, many system- 
atic uncertainties cancel. Thanks to these improvements in the experimen- 
tal method, the details and precision of the AIRFLY measurements surpass 
that of previous experiments. The application of the AIRFLY results to ultra 
high energy cosmic ray experiments which employ the fluorescence detection 
technique was also studied. A parameterization of the fluorescence yield as 
a function of altitude in terms of the measured relative band intensities and 
colhsional quenching pressures was derived. We estimated that the system- 
atic uncertainty on the cosmic ray shower energy associated with the pressure 
dependence of the fluorescence spectrum is reduced to a level of 1% by the 
AIRFLY results presented in this paper. 



9 Acknowledgements 



The authors wish to thank M. lanniUi, D. Pecchi, L. Sdogati, E. Tusi and 
G. Vitali of Universita di Roma Tor Vergata for the assistance in the con- 
struction of the pressure chamber, Ilya A. Shkrob and R.H. Lowers of the 
Argonne Chemistry Division for their support and assistance with the Van de 
Graaff maintenance and operations. This work was supported by the follow- 
ing agencies and organizations: United States Department of Energy, Office of 
Science, Offices of High Energy Physics and Basic Energy Sciences, under con- 
tract DE-AC02-06CH11357; BMBF grant no. 05 CU5VK1/8 and DFG grant 
no. KE 1151/1; project GA CR no. 202/05/2470 and project MSMT INGO; 
Istituto Nazionale di Fisica Nucleare, Italy; Center of excellence CETEMPS 
of Universita de I'Aquila. 



33 



References 



[I] R.M. Baltrusaitus et al, Nucl. Instrum. Meth. Phys. Res. A 240 (1985) 410. 
[2] T. Abu-Zayyad et al, Nucl Instrum. Meth. Phys. Res. A 450 (2000) 253. 

[3] J. Abraham et al, Nucl Instrum. Meth. Phys. Res. A 523 (2004) 50. 

[4] M. Fukushima, Prog. Theor. Phys. Suppl. 151 (2003) 206. 

[5] L. Scarsi et al, Proc. 27th Int. Cosmic Ray Conf., Hamburg (2001) 839 

[6] F. W. Stecker et al, Nucl. Phys. B136 (2004) 433 

[7] A. N. Bunner, Cosmic Ray Detection by Atmospheric Fluorescence, PhD thesis, 
Graduate School of Cornell University, February 1967. 

[8] G. Davidson and R. O'Neil, J. Chem. Phys. 41 (1964) 3946. 

[9] F. Kakimoto et al, Nucl. Instrum. Meth. Phys. Res. A 372 (1996) 527. 

[10] M. Nagano et al, Astropart. Phys. 20 (2003), 293; 22 (2004), 235. 

[II] J. W. Belz et al, Astropart. Phys. 25 (2006) 129. 

[12] P. Colin et al, |astro-ph/0612111 0, to appear in Astropart. Phys. 

[13] V. de Souza, G. Medina-Tanco, and J. A. Ortiz, Astropart. Phys. 25 (2006) 84. 

[14] B. Keilhauer et al, Astropart. Phys. 25 (2006) 259. 

[15] F. Arqueros et al, Astropart. Phys. 26 (2006) 231. 

[16] A. Lofthus and P. H. Krupenie, J. Phys. Chem. Ref. Data, 6 (1977) 113. 

[17] K. B. Mitchell, J. Chem. Phys. 53 (1970) 1795 

[18] J. M. Calo and R. C. Axtmann, J. Chem. Phys. 54 (1971) 1332 

[19] High Voltage Engineering Corporation model KS Van de Graaff. 

[20] An Eimac Y646 cathode electron source was used from Communications and 
Power Industries, Eimac Division, 607 Hansen Way, Palo Alto CA 94304 USA. 

[21] M. Conde et al, Proceedings of Particle Accelerator Conference, PAC05, 
Knoxville, Tennessee, USA, May 16-20,2005, TPAE016. 

[22] Bergoz Instrumentation, Espace Allondon Quest, 156, rue du Mont Rond, 01630 
Saint Genis Pouilly, France, model ICT-088-070-20:l. 

[23] Hamamatsu Photonics K. K., 314-5, Shimokanzo, Toyooka-village, Iwata-gun, 
Shizuoka-ken, 438-0193, Japan. 

[24] Varian, Inc., 121 Hartweh Avenue, Lexington, MA 02421, model SHOIOOIUNIV. 



34 



[25] Ocrlikon Leybold Vacuum GmbH, Bonner Strasse 498, Cologne D-50968, 
Germany, model Ceravac CTR 90. 

[26] Spectra-Physics, 150 Long Beach Boulevard, Stratford, Ct 06615, USA. 

[27] Andor Technology PLC, 7 Millenium Way, Springvale. Business Park, Belfast, 
BT12 7AL, Northern Ireland, model DV420 - BU2 CCD camera. 

[28] National Institute of Standards and Technology, Gaithersburg, MD, 20899, USA 

[29] A. R. Filippelli, S. Chung, and C. C. Lin, Phys. Rev. A29 (1984) 1709. 

[30] D. J. Eckstrom, J. S. Dickinson, and R. M. Hill, J. Chem. Phys. 82 (1985) 1792. 

[31] F. R. Gilmore, R. R. Laher, and P. J. Espy, J. Phys. Chem. Ref. Data, 21 (1992) 
1005. 

[32] J. T. Pons, R. S. Schappe, and C. C. Lin, Phys. Rev. A53 (1996) 2239. 

[33] W. L. Borst and. M. Imami, J. Appl. Phys. 44 (1973) 1133. 

[34] F. Blanco, and F. Arqueros, Phys. Lett. A345 (2005) 355 

[35] R. M. Sternheimer, S. M. Seltzer, and M. J. Berger, Phys. Rev. B26 (1982), 
6067. 

[36] S. Agostinein et al, Nucl. Instrum. Meth. A506 (2003) 250. 

[37] J. Baro' et al, Nucl. Instrum. Meth. BlOO (1995) 31; F. Salvat, J. M. Fernandez- 
Varea, and J. Sempau, Workshop Proceedings, Barcelona, Spain, 4-7 July 2006, 
NEA#06222, ISBN: 92-64-02301-1. 

[38] J. Belz et al, Astropart. Phys. 25 (2006) 57. 

[39] U.S. Standard Atmosphere 1976, U.S. Government Printing Office (Washington, 
D.C., 1976). 



35 



